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We reconsider the problem of optimal time to sell a stock studied by Shiryaev, Xu and Zhou 1 1] using path 
integral methods. This method allows us to confirm the results obtained by these authors and extend them to a 
parameter region inaccessible to the method used in Q]]. We also obtain the full distribution of the time tm at 
which the maximum of the price is reached for arbitrary values of the drift. 
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^ INTRODUCTION 

C/^ In the preceeding paper, A. Shiryaev, Z. Xu and X.Y. Zhou UJ] ask about the optimal tune to sell a stock over a certain time 

interval [0, T], knowing that the price is a geometrical Brownian motion with a certain average return over the risk-free rate a~r 
and a certain volatility a. The answer to this question depends on the value of the adimensional parameter a — {a — r)/a'^. 
The method used by the authors allow them to prove that whenever a > 1/2 the optimal selling time r* is always at the end 
of the interval, r* = T, whereas t* = in the case a < 0. In financial words, "good" stocks with a sufficiently large average 
return should be sold as late as possible, whereas one should immediately get rid of "bad stocks". These results are clearly very 
interesting; however one feels unsatisfied by the fact that the authors' method do not allow them to treat the case < a < 1/2. 
Ch They discuss this point in the conclusion, mentioning (a) a working paper [2] based on an alternative method showing that one 
I should in fact sell immediately as soon as a < 1/2 and (b) that the case a < 1/2 is not interesting financially because "most 
stocks realize a > 1/2 by a large margin". 
I— I , The aim of this short note is to reconsider the problem using path integral methods which are well known in physics but 
^ ' perhaps less well known in financial mathematics. This method allows one to treat all values of a on the same footing. We 
confirm the results of 01 and extend them to the < a < 1/2 interval. In fact, we show that there is an exact symmetry 
QQ in the problem that relates the problem with a > 1/2 to the problem with a < 1/2. Our method furthermore allows us to 
f — ■ garner additional results, such as the distribution of the time at which the maximum of the price is reached. We find that this 
OO distribution has inverse square root singularities both at t„i — and tm = T for all values of a; however, the amplitude of the 
divergence at t„i = is stronger when a < 1/2 and weaker when a > 1/2. This gives a more precise picture to the results of 
0^ Shiryaev, Xu and Zhou. For a = 1/2, the problem is degenerate and the two peaks have exactly the same amplitude (in fact, the 
distribution is symmetric under T — 

Finally, we do not agree with the statement that a < 1/2 is not interesting financially. The numbers provided in fl'] are 
based on the S&P500 index returns, and are therefore much too optimistic: first, there is an obvious selection bias since badly 
performing stocks leave the index; second, the volatility of the index is two to three times smaller than the volatility of individual 
stocks, thanks to the diversification effect. An annualized volatility above 40% is in fact not uncommon, in particular for small 
;_3 to medium caps - whereas the S&P500 only includes large caps. Fig. 1 shows the time series of the (implied) S&P500 index 
Ci volatility and the average stock volatility, in the period 2000-2007. With an average annual return of 10%, an interest rate of 5% 
and a volatility of 40% annual, the parameter a is found to be 0.3125 < 1/2. 

We hope that this short note will shed a useful light on the work of A. Shiryaev, Z. Xu and X.Y. Zhou, and that it will convince 
the reader that path integral methods are extremely powerful to solve a variety of random walk problems. We refer the reader to 
a short review paper by one of us iJt] on this topic, see also 
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THE SET-UP 

In this section we give the set-up of the problem using physicists notations. We assume, as in IH, that the price Pt of a stock 
follows a geometric Brownian motion: 

Pt = exp[^i + crBt], (1) 

where /i is Ito corrected drift and Bt the standard Brownian motion. We will use below the notation x{t) = fit + aBt for the 
drifted Brownian motion, with by convention x{Q) = 0. In Ref. yj, the authors introduce the notation a — iijo^ + 1/2. A 
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FIG. 1: Time series of the (implied) S&P500 index volatility (the VIX, bottom green curve) and the average implied stock volatility for SPX 
stocks (middle red curve) and mid-cap stocks (upper black curve), in the period 2000-2007. It is clear that stock volatilities are on average two 
to three times larger than the volatility of the index, and that values of a above 40% are not uncommon. 

'good' stock in the financial language corresponds to a positive drift /i > (i.e., a > 1/2) and a 'bad' stock coiTesponds to a 
negative drift /i < (a < 1/2). In terms of the real return a of the stock, the condition /i > translates into a — r > (t^/2, 
where a — r is the excess return over the risk-free rate r. Note that the process that we talk about is the real world process and 
not the risk-neutral one, which has no meaning for the question raised in Ref. [1]. 

Let us consider the evolution of the stock price over a fixed time interval < i < T. It is intuitively obvious that the maximum 
of a drifted Brownian motion and hence that of the stock price Pt is most likely to occur att = T (for /i > 0) and t ^ (for 
/i < 0). Thus, it obviously makes sense to sell a 'good' stock (/i > 0) at the end of the interval t = T, whereas a 'bad' stock 
(/i < 0) at the begining of the interval t = 0. This intuitive results are put on a more rigorous mathematical footing in the rest 
of this note by (i) calculating exactly, using path integral methods, the maximal relative error as defined in Ref. f^, but for all 
values of /i and (ii) also by computing the full probability density of the time tm at which the maximum of the price occurs for 
all fx. 

Let Mx denote the maximum price of the stock over the interval [0, T], i.e., 

Mt = max Pt. (2) 

0<t<T 

Evidently, the optimal time to sell the stock is the one where the difference between the price of the stock and its maximal 
value Mt is minimal. A convenient way to estimate this optimal time is to consider the relative eiTor at a fixed time r where 

< T < T 

-<-)-(^)— (^) 

where E denotes the expectation value over all realizations of the Brownian motion. Minimizing (t, T) over < r < T gives 
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the optimal time r* . In other words, r* is the time at which the ratio 

5^(T,T) = l-,v(r,r) = E(^^) (4) 
is maximal. The goal is to estimate S^{t, T) and then maximize it with respect to < t < T. Using the trivial identity 



Mt = max Pr — max [exp(a;(T))l — exp 

0<T<T 0<T<T 

one can rewrite (r, T) in Eq. (|4]i as 

S,,{t,T) = E [exp (-{Mt - x{t)} 

where 



max x(t) 

0<T<T 



(5) 



(6) 



Mt = max x(t) = InfMrl (7) 

is the maximum of the drifted Brownian motion x{t) over < t < T. Note that throughout this paper, we will use t as the 
running time and t as a fixed time. 

Let us consider the random variable ?/(r) = Mt — x{t) at a fixed time r and let Pfi{y, r) denote its probability density 
function (pdf). Once we know r), then from Eq. Q, we can evaluate 



S^{t,T)= / dye-y P^iy,T). (8) 
Jo 

To evaluate the pdf (y, r), we need the joint pdf of Mt and x{t) at fixed r. 

JOINT DISTRIBUTION OF AIt AND a:(r) 

It is convenient to compute first the cumulative probability 

Ff,{x, TO, t) = Prob[a:(T) = x and Mt < m] (9) 

where the walk starts at the origin x{Q) ~ and Mt is the global maximum of the walk in [0, T]. This cumulative probability 
can be computed using a path-integral approach as detailed below. 

Clearly F^{x, m, r) is the probability that a drifted Brownian motion x{t) in < t < T, starting from a;(0) = 0, reaches 
x{t) = a; at a fixed time t — t and in addition, stays below the level m for all < t < T. The last condition comes from the 
fact that if the global maximum Mt < to, the path necessarily stays below the level to for all < i < T. An example of such 
a path is seen in Fig. 2. To compute F^{x, to, t), it is convenient to consider the shifted process y{t) = rn — x{t) so that the 
process y{t) evolves, as 

dy — —dx = —fidt — adBf. (10) 

Thus the shifted process y{t) represents a Brownian motion with a drift — /i, opposite to that of x{t). In terms of the process 
y{t), F^{x, to, r) is just the probability that the process y{t), starting at y(0) — to, reaches the point j/(t) — m — xatt — t and 
stays positive in the whole interval < i < T. An example of such an event is shown in Fig. 3 

For the process y(t) in Eq. dTol l, let us first define the propagator G^^{y, yo, t) that denotes the probability that the process 
starting at yo?itt ^ 0, reaches y at time t, but staying positive in between, i.e., in [0, t]. One can then easily express F^{x, to, t) 
in terms of this propagator as (see Fig. 3) 

F^{x,m,T)=Gt^{m-x,m,T) / G+{y',m~x,T-T)dy'. (11) 

Jo 

In writing Eq. (fTTl l. we have spht the interval [0, T] into two parts: [0, r] and [r, T]. In the first interval (see Fig. 3), the process 
propagates from the initial position y{0) = to to ^(t) = to — x in time r (staying positive in between), hence explaining the first 
factor G^^{m — x, to, r) in Eq. (fTTT i. In the second interval, the process starting at the new 'initial' position to — x, propagates 
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FIG. 2: A realization of the drifted Brownian motion x{t) in t G [0, T], starting at x{0) — 0, reaching x{t) = x ait — t and staying below 
the level m for all < t < T. 



to a final position y' in time T — t, staying positive in between. Also, the final position y' can be any positive number and one 
has to integrate over it. This explains the second factor in Eq. (fTTl i. Of course, in writing the path decomposition in Eq. (fTTT l we 
have used the renewal property of a Brownian motion (valid due to its Markovian nature) which implies that the two intervals 
(left of T and right of r) are statistically independent. 

Evaluation of the propagator G^^{y, yo,T): Using a physicist interpretation of Eq. ( fTOl i. we note that the Langevin noise 
ri{t) — dBt/dt is a Gaussian white noise with the associated measure, Prob [{?7(/:)}, < t < t] (x exp [— ^ r]^{t)dt\. 
Substituting, 77 (t) ~ {y + ^i)/a from Eq. ( fTOb . one can express the propagator as a path integral 



G-a{y,yo,T) = 



v{T)=y 



y{o)=yo 



Vy{t) exp 



1 

2^ 



dt (y + fiY 



.t=0 



(12) 



where [nt=o ^iui^))] indicator function that enforces the path to stay positive in the interval t G [0, r]. The rhs of Eq. (fT2l l 
can be rearranged (by expanding the square {y + fi)'^ and performing the time integral) as 



exp 



2a2 



:{y-yo) 



Go{y,yo,T) 



(13) 



where G^(?/, yo, t) is the propagator associated with the driftless (/i = 0) Brownian motion 



Go{y,yo, 



v{t)=V 



y{o)=yo 



Vyit) exp 



1 

'2^ 



dtf 



U^^iyit)) 



(14) 
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FIG. 3: A realization of the shifted Brownian motion y{t) with drift —fj, in t G [0, T], starting at y{0) = m, reaching y(r) = m — xatt = T 
and staying positive for all < t < T. 



This propagator, which denotes the probabiUty that a driftless Brownian motion propagates from yo to y in time r without 
crossing the origin in between, can be evaluated very simply by the standard method of images Isll^ or alternatively by the path 
integral method [3J and has the well known expression 



G(j"(y,2/o,T) 



1 



exp 



jy - yo) 



exp 



jy + yof 

2<7^T 



(15) 



Substituting this in Eq. ( fTSl l. one then has the required propagator. Using this explicit expression for Gl^ (y, yo, r) one can also 
easily evaluate the following integral 



1 

Glf,{y,yo,T)dy = - 



erfc 



yo - 



r 



exp 



'2^J.yo 



erfc 



yo + fJ-T 



(16) 



where erfc(a;) = e " du is the complementary error function. Assembling these results in Eq. (fTTI) gives us an explicit 

expression for the cumulative probability 



Ffj,{x,m,T) 



e 2<7^ <T^ 



2V2TTa^T 



^2 (2to---x)^ 
e 2<t2x — e 2<t2x 



erfc 



^(T — t) \ 2m (m - a;) 



erfc 



m ^ X + /i(T — r) 



(17) 

The joint pdf Q^{x, m) of x(t) = a; at fixed r and Af^ = ^ can then be obtained by taking the derivative of to, t) with 
respect to to, i.e.. 



,(a;(T) = x, Mt = to) 



dF^{x,m,T) 
dm 



(18) 
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EVALUATION OF THE RELATIVE ERROR r^(r, T) 



Having obtained the joint pdf Qfj,{x{T) — x, Mt — m) in Eqs. (fTSl l and (fTTb . we can easily find the pdf -P^(y, r) of the 
variable y = Mt — x 



Qfj,{x, m) S {y — (m — x)) dx dm 
Qp(m - y,m) dm. 



The above integral can be performed exactly (we skip the details here). One obtains the following expression 



1 



where 



9f^{y,T) 



1 



exp 
erfc 



(tV 2 



jy + IJ-T? 



VV2^ 



V\/2^ 



Note from the explicit expression of r) the following symmetry 

P^(2/,r)=P_^(y,r-T) 



(19) 
(20) 

(21) 
(22) 

(23) 



which has the simple physical meaning of time-reversal symmetry, i.e., when the process propagates in the reverse time direction, 
one gets the same measure provided one also reverses the sign of the drift /x. 

Having obtained the pdf Pf_i{y, t), one can then evaluate the relative error r^(r, T) = 1 — S'^ (r, T) where 



/•oo 

5^(r,r)= / dye-yp^{y,T) 
Jo 



Evidently S'^ (r, T) also has the same time-reversal symmetry namely 

S^iT,T)^S^^{T~T,T) 



(24) 



(25) 



While it is difficult to do the integral in Eq. ( |24] | analytically, one can easily evaluate it using Mathematica. Besides, the general 
feature of (r, T) as a function of r can be inferred by just studying the asymptotic properties of the integral in Eq. (l24b in the 
limit T ^ and r ^ T. In Fig. 4, we show a plot of S^{t, T) for three different values of ^ = 0.1, /i = and /i = —0.1 upon 
setting T = 1 and a = 1. 

Optimal Time r* : To find the optimal time r* we need to minimize (r, T), i.e., maximize 5^ (r, T) with respect to r e [0, T] . 
It is evident from Fig. 4 and also from the expression of ^^(t, T) that for all values of /i, S^{t, T) has two local maxima at the 
endpoints of the interval [0, T], i.e., respectively at t = and t — T. However, for fi > 0, the maximum at r = T has a larger 
value implying that for fi > 0, t* = T. By the symmetry manifest in 5;^(r, T) — S'_p(r — r, T) it follows that for /i < 0, the 
maximum at r = has a higher value implying r* = for /i < 0. Exactly at /i = 0, both local maxima at r = and t — T 
have the same value (S'o(t, T) is completely symmetric around the midpoint t = T/2) implying that for /i = 0, both r* = 
and T* — T are optimal. 

The optimal value S'^ (r* , T) is actually easier to evaluate since for t = or r = T (at the end-points), the integral in Eq. 
(|24] | can be carried out explicitly. Omitting details of this integration, we get the following expression for the optimal relative 
error for all /i 



r(r*,r) - 1-S^{t*,T) 




2|m| 



exp 




(26) 



Note that the optimal relative error is evidently a symmetric function of /i as manifest in the above result. 
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FTG. 4: Plots of S';j(t, T) vs. r obtained from Eq. ( 124b for three different values of the drift jj. 

r = 1 and o- = 1. The symmetry S^{t, T) = ^-^(r - r, T) is evident. 



0.1, fi — and — —0.1. We have set 



In the preceeding paper |1], Shiryaev, Xu and Zhou also obtained an expression of the optimal relative error r(T*,T) by a 
completely different method. Their notations are slightly different from above. In their notation, yu = (a — 1 /2)o'^ and also their 
result for r{T* ,T) is in terms of the probability distribution of a Gaussian random variable with zero mean and unit variance, 
= -h= e^" /^du which is related to the complimentary eiTor function via 



■ erfc 



V2) 



(27) 



However, their method allows them to obtain an explicit expression for the optimal relative error only in the range a > 1/2 (i.e., 
/i > 0) and a < (i.e., /i < — a^/2). In these ranges, their expressions for the optimal relative error (Eqs. (9) and (1 1) in fill) 
reduce precisely to our compact result in Eq. ( |26] l, upon identifying /i = (a — l/2)cr^ and $(a;) as in Eq. (|27T l. However, they 
do not have any result in the range < a < 1/2, i.e., for — cr^/2 < < 0. In contrast, our result in Eq. ( |26] l is valid for all fi 
(and hence for all a) and is therefore more general. In addition, their method somehow does not detect the symmetry of r (t* , T) 
under fi —fi which is manifest in our path integral approach. 
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THE EXACT DISTRIBUTION OF THE TIME OF THE OCCURRENCE OF THE MAXIMUM FOR A BROWNIAN MOTION 

WITH DRIFT fi 



Minimizing the relative error r(T, T) with respect to t is one way of estimating the optimal time r* at which one should sell 
a stock over a fixed investment time horizon T, as explained above. Another alternative and direct measure would be to first 
derive the probability density p(tm, T) of the time at which the maximum Mt of a stock price over [0, T] actually occurs. 
This density p{tm, T) will typically have a peak (or more peaks). The value of i,„ = t* at which the strongest peak of p{t„i, T) 
occurs can then be taken as an alternative measure for the optimal time to sell a stock, since the maximum of the price is most 
likely to occur at tm = t* . 

In this section, we compute exactly the density p^{t„i, t) of tm for a Brownian motion x{t) with drift /i. Since the stock price 
Pt — exp[x{t)] is just the exponential of x{t) under the Black-Scholes scenario, the maximum Mt of the stock price Pt occurs 
exactly at the same time <,„ where x{t) itself achieves its maximum. For the case fi — 0, the density po{tm, T) was computed 
by Levy and is given by the derivative of an arcsine form, i.e.. 



(28) 



Recently, using an appropriate path integral method, the density of tm was computed exactly for a Brownian motion up to its 
first-passage time |8J and also for a variety of constrained Brownian motions such as Brownian excursions, Brownian bridges, 
Brownian meanders etc. 101 ■ Here we adapt this path integral method to compute the density p^{tm, T) for a Brownian motion 
with arbitrary drift /i. 

To compute the density Pfj,{tm,T) the strategy would be to first compute the joint density of tm as well as the maximum 
Mt — m itself, i.e., p^{tm, m, T) and then integrate over ni to obtain the marginal density, p^{tm, T) — Pf^{tm,'m, T) dm. 
The joint density p^{tm, rn, t) is the proportional to the sum of the statistical weights of all paths that start at the origin x(0) = 0, 
reaches the value x{tm) — rn for the first-time att — tm and then stays below the level m at all subsequent times up to r,i.e., in 
the interval [tm, T]. To enforce the conditions that x{t) < m in the two intervals t e [0, tm] and t £ [tm, T] and that exactly at 
tm the path reaches x{tm) = rn, poses a problem for a continuous-time Brownian motion. This is because if a Brownian motion 
crosses a level m at a given time tm then it must cross and re-cross the same level m an infinite number of times in the vicinity 
of t = tm- Hence it is impossible to enforce the above constraints simultaneously for a continuous-time Brownian motion. Note 
that for lattice random walks this does notpose any problem. To get around this difficulty with the continuous-time Brownian 
motion, one introduces a small cut-off e ||4|,|8|], i.e., one assumes that the path, starting at x(0) — reaches the level to — e at 
time tm, staying below to for all < i < and then starting at m — e at t = tm stays below the level m for all tm < t < T 
(see Fig. 5 for such a realization). Finally one takes the limit e ^ at the end of the calculation. 

Comparing Figs. (2) and (5), it is clear that the paths that contribute to the joint probability density Pfj^{tm,'m,T[e) are 
identical to those that contribute to F^{x, to, r) with the replacements x = m — e and r = i„, in Eq. ( [TtI i. i.e., p^{tm,'ni, T|e) oc 
F^{m — e, TO, tm)- Substituting x — ra — e and t = tm in Eq. ( [TtI i and taking the e ^ limit we find, to leading order in e. 



Pf_,{tm,m,T[e) 



y/2TTa^{T~tm) 




T-tr. 



(29) 



where the constant of proportionaUty A, which is function of e, is determined from the normalization, 
/q dtm dmp^(tm,'m,T[e 0) = 1. This fixes A = a^/e^. Integrating T) (now the cut-off e has been 

set to 0) over to finally gives the marginal density p^ {tm , T) in a closed form 



Pp,(tm,T) = 



1 



Tr\/tmiT — tm) 



h(tm,tJ.)h{T -tm,-tJ-) 



(30) 



where 



h{tm,lj) = exp - 



■ntrr 



erfc 




(31) 



The density p^{tm, T) given in Eqs. (l30b and (ISTT l is the main result of this section. Evidently, for /i = 0, one recovers from 
this the well known arcsine result of Levy in Eq. ( |28] |. Note that the density p^{tm, T) also has a symmetry similar to that in 
Eq. (|25] | namely 



Pf^itm,T) ^p^i_,{T -tm,T). 



(32) 
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FIG. 5: A realization of the drifted Brownian motion x{t) in t G [0, T], starting at 2;(0) = 0, reaching x{tm) — m — e at t — tm and staying 
below the level m for all < t < T. 



This symmetry is also evident in Fig. 6 where we plot the density p^{t,n, T) in Eq. (|30] l for ^ — 1, ji — Q and /i = — 1 upon 
setting T = 1 and a = 1. 

We note from Eq. ( |30] l as well as from Fig. 6 that for all values of /i, the density p^{tm, T) has two peaks (actually has square 
root divergences) at the two end points tm — ^ and tm = T, 



where the amplitude 



exp 



»0,T) 

T,r) 

2a2 



A ,{T) 



(33) 
(34) 

(35) 



However, for > 0, the divergence at — T is stronger than that at t„i — since y4_p(T) > ^^(T). On the other hand, 
for /i < 0, the opposite is true. At /i = 0, both ends have the same divergences as the density is completely symmetric around 
t„i — T/2. Thus, we conclude that the maximum of the Brownian motion with drift fi is most likely to occur at t„i = T for 
/J, > 0, at t„i = for < 0, and for /i = both t„i = and i,„ ~ T are equally likely. This then leads us to identify the optimal 
time T* to sell the stock (within the black-Scholes economy model) to be r* = T for /i > 0, r* = for /i < 0, and t* = 0, T 
(equally likely) for fi — Q. Thus, based on the analysis of the density {tm , T) we draw the same conclusion as was obtained 
from the optimization of the relative error in the previous sections. 
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FIG. 6: Plots of Pfj.{tm, T) vs. tm for three different values of the drift fi — 1.0, fi = and jj, = —1.0. We have set T = 1 and a — 1. The 
symmetry p^{tm,T) = P~^{T ~ tm,T) is evident. 
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